Simplifying Genotyping of Mutants from Genome Editing with a Parallel qPCR-Based iGenotype Index.

Targeted genome editing is a powerful tool in reverse genetic studies of gene function in many aspects of biological and pathological processes. The CRISPR/Cas system or engineered endonucleases such as ZFNs and TALENs are the most widely used genome editing tools that are introduced into cells or fertilized eggs to generate double-strand DNA breaks within the targeted region, triggering cellular DNA repair through either homologous recombination or non-homologous end joining (NHEJ). DNA repair through the NHEJ mechanism is usually error-prone, leading to point mutations or indels (insertions and deletions) within the targeted region. Some of the mutations in embryos are germline transmissible, thus providing an effective way to generate model organisms with targeted gene mutations. However, point mutations and short indels are difficult to be effectively genotyped, often requiring time-consuming and costly DNA sequencing to obtain reliable results. Here, we developed a parallel qPCR assay in combination with an iGenotype index to allow simple and reliable genotyping. The genotype-associated iGenotype indexes converged to three simple genotype-specific constant values (1, 0, -1) regardless of allele-specific primers used in the parallel qPCR assays or gene mutations at wide ranges of PCR template concentrations, thus resulting in clear genotype-specific cutoffs, established through statistical analysis, for genotype identification. While we established such a genotyping assay in the Xenopus tropicalis model, the approach should be applicable to genotyping of any organism or cells and can be potentially used for large-scale, automated genotyping.

Muscle Fiber Characteristics and Transcriptome Analysis in Slow- and Fast-Growing Megalobrama amblycephala.

Growth is an important trait in aquaculture that is influenced by various factors, among which genetic regulation plays a crucial role. Megalobrama amblycephala, one of the most important freshwater species in China, exhibits wide variations in body mass among individuals of the same age within the same pool. But the molecular mechanisms underlying wide variation in body mass remain unclear. Here, we performed muscle histological and transcriptome analysis of muscle tissues from Fast-Growing (FG) and Slow-Growing (SG) M. amblycephala at the age of 4 months old (4 mo) and 10 months old (10 mo) to elucidate its muscle development and growth mechanism. The muscle histological analysis showed smaller diameter and higher total number of muscle fibers in FG compared to SG at 4 mo, while larger diameter and total number of muscle fibers were detected in FG at 10 mo. The transcriptome analysis of muscle tissue detected 1171 differentially expressed genes (DEGs) between FG and SG at 4 mo, and 718 DEGs between FG and SG at 10 mo. Furthermore, 44 DEGs were consistently up-regulated in FG at both 4 mo and 10 mo. Up-regulated DEGs in FG at 4 mo were mainly enriched in the pathways related to cell proliferation, while down-regulated DEGs were significantly enriched in cell fusion and muscle contraction. Up-regulated DEGs in FG at 10 mo were mainly enriched in the pathways related to cell proliferation and protein synthesis. Therefore, these results provide novel insights into the molecular mechanism of M. amblycephala muscle growth at different stages, and will be of great guiding significance to promote the fast growth of M. amblycephala.

Predictive value of immunotherapy-induced inflammation indexes: dynamic changes in patients with nasopharyngeal carcinoma receiving immune checkpoint inhibitors.

BACKGROUND: Immune checkpoint inhibitors (ICIs) have achieved substantial advancements in clinical care. However, there is no strong evidence for identified biomarkers of ICIs in NPC. METHODS: In this retrospective study, 284 patients were enrolled into a training or validation cohort. Inflammatory indexes based on peripheral blood parameters were evaluated, including the systemic immune-inflammation index (SII), the neutrophil-lymphocyte ratio (NLR), the platelet-lymphocyte ratio (PLR), the lymphocyte-to-C-reactive protein ratio (LCR), and the lymphocyte-monocyte ratio (LMR). The optimum cut-off value for patient stratification was identified using X-tile. The Kaplan-Meier method and Cox's proportional regression analyses were used to identify prognostic factors. RESULTS: Immunotherapy significantly changed the levels of SII, NLR, PLR, LCR and LMR in NPC patients. Patients with lower SII, NLR, and PLR, as well as those with higher LCR and LMR, before immunotherapy had superior PFS (all p < 0.05). Moreover, PFS in the decreased SII, reduced NLR and increased LMR group was significantly longer than in the opposite group (all p < 0.05). Both univariate and multivariate analyses validated that baseline SII and LMR, and the immunotherapy-related SII reduction and LMR elevation were independent prognostic factors for PFS in advanced NPC patients receiving ICIs. CONCLUSIONS: Immune checkpoint inhibitor treatments significantly changed the levels of SII, NLR, PLR, LCR and LMR in NPC patients treated with immunotherapy. A lower baseline SII and a higher baseline LMR, and a reduction in SII and an elevation in LMR after immunotherapy are favorable factors for predicting survival among advanced NPC patients.

There is no strong evidence for identified biomarkers of immune checkpoint inhibitors (ICIs) in nasopharyngeal carcinoma (NPC).Lower baseline SII and higher baseline LMR were related to better PFS. The dynamic changes of SII and LMR were independent prognostic factors for the survival of NPC patients receiving ICIs.Neutrophils, platelets, lymphocytes, and monocytes can be used as cheap and valuable biomarkers for predicting tumor response in NPC on immunotherapy.

HSF1 protects cells from cadmium toxicity by governing proteome integrity.

BACKGROUND: Cadmium toxicity has been associated with disruption of protein homeostasis by interfering with protein folding processes. Heat shock factor 1 (HSF1) coordinates the rapid and extensive cellular response to maintain proteomic balance facing the challenges from many environmental stressors. Thus, we suspect that HSF1 may shield cells from cadmium toxicity by conserving proteome integrity. RESULTS: Here, we demonstrate that cadmium, a highly poisonous metal, induces aggregation of cytosolic proteins in human cells, which disrupts protein homeostasis and activates HSF1. Cadmium exposure increases HSF1's phosphorylation, nuclear translocation and DNA bindings. Aside from this, HSF1 goes through liquid-liquid phase separation to form small nuclear condensates upon cadmium exposure. A specific regulatory domain of HSF1 is critical for HSF1's phase separation capability. Most importantly, human cells with impaired HSF1 are sensitized to cadmium, however, cells with overexpressed HSF1 are protected from cadmium toxicity. Overexpression of HSF1 in human cells reduces protein aggregates, amyloid fibrils and DNA damages to antagonize cadmium toxicity. CONCLUSIONS: HSF1 protects cells from cadmium toxicity by governing the integrity of both proteome and genome. Similar mechanisms may enable HSF1 to alleviate cellular toxicity caused by other heavy metals. HSF1's role in cadmium exposure may provide important insights into the toxic effects of heavy metals on human cells and body organs, allowing us to better manage heavy metal poisoning.

Identification of thyroid hormone response genes in the remodeling of dorsal muscle during Microhyla fissipes metamorphosis.

Introduction: Extensive morphological, biochemical, and cellular changes occur during anuran metamorphosis, which is triggered by a single hormone, thyroid hormone (TH). The function of TH is mainly mediated through thyroid receptor (TR) by binding to the specific thyroid response elements (TREs) of direct response genes, in turn regulating the downstream genes in the cascade. The remodeling of dorsal skeletal muscle during anuran metamorphosis provides the perfect model to identify the immediate early and direct response genes that are important during apoptosis, proliferation, and differentiation of the muscle. Methods: In our current study, we performed Illumina sequencing combined with single-molecule real-time (SMRT) sequencing in the dorsal muscle of Microhyla fissipes after TH, cycloheximide (CHX), and TH_CHX treatment. Results and Discussion: We first identified 1,245 differentially expressed transcripts (DETs) after TH exposure, many of which were involved in DNA replication, protein processing in the endoplasmic reticulum, cell cycle, apoptosis, p53 signaling pathway, and protein digestion and absorption. In the comparison of the TH group vs. control group and TH_CHX group vs. CHX group overlapping gene, 39 upregulated and 6 downregulated genes were identified as the TH directly induced genes. Further analysis indicated that AGGTCAnnTnAGGTCA is the optimal target sequence of target genes for TR/RXR heterodimers in M. fissipes. Future investigations on the function and regulation of these genes and pathways should help to reveal the mechanisms governing amphibian dorsal muscle remodeling. These full-length and high-quality transcriptomes in this study also provide an important foundation for future studies in M. fissipes metamorphosis.

Transcriptome profiling reveals gene regulation programs underlying tail development in the Ornamented Pygmy frog Microhyla fissipes.

Introduction: Tadpole tail develops from the tailbud, an apparently homogenous mass of cells at the posterior of the embryo. While much progress has been made in understanding the origin and the induction of the tailbud, the subsequent outgrowth and differentiation have received much less attention, particularly with regard to global gene expression changes. Methods: By using RNA-seq with SMRT and further analyses, we report the transcriptome profiles at four key stages of tail development, from a small tailbud to the onset of feeding (S18, S19, S21 and S28) in Microhyla fissipes, an anuran with a number of advantages for developmental and genetic studies. Results: We obtained 48,826 transcripts and discovered 8807 differentially expressed transcripts (DETs, q < 0.05) among these four developmental stages. We functionally classified these DETs by using GO and KEGG analyses and revealed 110 significantly enriched GO categories and 6 highly enriched KEGG pathways (Protein digestion and absorption; ECM-receptor interaction; Pyruvate metabolism; Fatty acid degradation; Valine, leucine and isoleucine degradation; and Glyoxylate and dicarboxylate metabolism) that are likely critically involved in developmental changes in the tail. In addition, analyses of DETs between any two individual stages demonstrated the involvement of distinct biological pathways/GO terms at different stages of tail development. Furthermore, the most dramatic changes in gene expression profile are those between S28 and any of the other three stages. The upregulated DETs at S28 are highly enriched in "myosin complex" and "potassium channel activity", which are important for muscle contraction, a critical function of the tail that the animal needs by the end of embryogenesis. Additionally, many DETs and enriched pathways discovered here during tail development, such as HDAC1, Hes1 and Hippo signaling pathway, have also been reported to be vital for the tissue/organ regeneration, suggesting conserved functions between development and regeneration. Conclusion: The present staudy provides a golbal overview of gene expression patterns and new insights into the mechanism involved in anuran tail development and regeneration.

TIP30 overcomes gefitinib resistance by regulating cytoplasmic and nuclear EGFR signaling in non-small-cell lung cancer.

Epidermal growth factor receptor (EGFR) tyrosine kinase inhibitors (TKIs) (eg, gefitinib) exert potent therapeutic efficacy in non-small-cell lung cancer (NSCLC) harboring EGFR-activating mutations. However, the resistance to EGFR TKIs limits their clinical therapeutic efficacy. TIP30, a newly identified tumor suppressor, appears to be involved in the regulation of cytoplasmic and nuclear EGFR signaling in NSCLC. Our previous study demonstrated that TIP30 regulated EGF-dependent cyclin D1 transcription in human lung adenocarcinoma and suppressed tumorigenesis. In the present study, the involvement of TIP30 in combating gefitinib resistance in NSCLC was determined for the first time in vitro and in vivo. Gain and loss of function studies showed that overexpression of TIP30 effectively sensitized cells to gefitinib in vitro, whereas TIP30 inhibition promoted gefitinib cell resistance. Moreover, TIP30 negatively regulated the activation of the p-AKT and p-MEK signaling pathways in PC9/GR. Importantly, PC9/GR harbored high levels of nuclear EGFR, and overexpression of TIP30 restored irregular EGFR trafficking and degradation from early endosomes to the late endosomes, decreasing the nuclear accumulation of EGFR, which may partly or totally inhibit EGFR-mediated induction of c-Myc transcription. Xenographic tumors induced by overexpression of TIP30 by PC9/GR cells in nude mice were suppressed compared with their original counterparts. Overall, it was revealed that TIP30 overexpression restored gefitinib sensitivity in NSCLC cells and attenuated the cytoplasmic and nuclear EGFR signaling pathways and may be a promising biomarker in gefitinib resistance in NSCLC.

The Changes in the Frog Gut Microbiome and Its Putative Oxygen-Related Phenotypes Accompanying the Development of Gastrointestinal Complexity and Dietary Shift.

There are many examples of symbiotic and reciprocal relationships in ecological systems; animal gut microbiome-host interactions are one such kind of bidirectional and complex relationship. Here, we utilized several approaches (16S rRNA gene sequencing, metagenomics, and transcriptomics) to explore potential gut microbiome-host interactions accompanying the development of gastrointestinal complexity and a dietary shift from metamorphosis to maturity in ornamented pygmy frogs (Microhyla fissipes). We identified the possible coevolution between a particular gut microbial group (increased putative fat-digesting Erysipelotrichaceae and chitin-digesting Bacteroides and Ruminococcaceae) and the host dietary shift [from herbivore to insectivore (high proportion of dietary chitin and fat)] during metamorphosis. We also found that the remodeling and complexity of the gastrointestinal system during metamorphosis might have a profound effect on the gut microbial community (decreasing facultative anaerobic Proteobacteria and increasing anaerobic Firmicutes) and its putative oxygen-related phenotypes. Moreover, a high proportion of chitin-digesting bacteria and increased carbohydrate metabolism by gut microbiomes at the climax of metamorphosis would help the frog's nutrition and energy needs during metamorphosis and development. Considering the increased expression of particular host genes (e.g., chitinase) in juvenile frogs, we speculate that host plays an important role in amphibian metamorphosis, and their symbiotic gut microbiome may help in this process by providing the nutrition and energy needs. We provide this basic information for the amphibian conservation and managements.

Metabolic switch in energy metabolism mediates the sublethal effects induced by glyphosate-based herbicide on tadpoles of a farmland frog Microhyla fissipes.

Glyphosate-based herbicides (GBHs) are widely-used agricultural chemicals, bringing potential detriments to aquatic organisms. Currently, our understanding of sublethal effects and underlying toxicologic mechanisms of GBHs are still limited, especially in amphibians. Here, the sublethal effects of a commercial GBH (KISSUN®) on tadpoles of a farmland dwelling frog, Microhyla fissipes, were investigated. The 10-d LC50 of "KISSUN®" GBH was 77.5 mg/L. Tadpoles exposed to 60-120 mg/L showed increased preference to higher temperature. After 10 days exposure, obvious growth suppression was observed in survived GBH-stressed tadpoles, characterized by dosage depended decrement in body mass, body width, total length, etc. GBH-stressed tadpoles also showed decreased tail length/snout-vent length ratio and smaller tail muscle fiber diameter. Comparative transcriptomics (control, 60 mg/L and 90 mg/L groups) was conducted to analyze the underlying molecular processes. GBH-stressed tadpoles showed downregulated transcription of ribosomal proteins and cytoskeleton proteins, which could explain their suppressed whole body and tail muscle growth. Moreover, GBH-stressed tadpoles showed transcriptional downregulation of carbohydrate and lipid catabolism, but upregulation of amino acid catabolism. It suggested a metabolic switch from carbohydrate and lipid to amino acid in these tadpoles. Accordingly, there was a trade-off between protein synthesis and energy production in respect to amino acid allocation, and it provided a metabolic explanation for why protein synthesis was downregulated and growth was suppressed in GBH-stressed tadpoles. In combination with existing literatures, we speculated that GBH might directly target the enzymes in carbohydrate and lipid catabolism, and this metabolic effect of GBH might be common to fish and amphibians. In conclusion, our study provided a systematic insight into the sublethal symptoms of GBH-stressed tadpoles, and a metabolic switch from carbohydrate and lipid to amino acid likely underlay some common toxic symptoms of GBHs on both fish and tadpoles.

MDM2 inhibitor RG7388 potently inhibits tumors by activating p53 pathway in nasopharyngeal carcinoma.

Nasopharyngeal carcinoma (NPC) is a high-risk head and neck cancer with poor clinical outcomes and insufficient treatments. The mouse double minute 2 homolog (MDM2) is the main molecular target in the clinical treatment of cancer. Indeed, MDM2 negatively regulates p53 through ubiquitin-dependent degradation. Thus, inhibition of MDM2-p53 interaction is a potential strategy for treating NPC. The latest generation MDM2 inhibitor, RG7388, shows increased potency and improved bioavailability compared to previous treatments. In this study, we investigated the efficacy and specificity of this inhibitor in NPC cell lines, and tumor-bearing mice were used to examine the therapeutic efficacy and effects of RG7388 treatment. The results showed that RG7388 potently decreased cell proliferation and activated p53-dependent pathway, resulting in cell cycle arrest and apoptosis. RG7388 significantly inhibited tumors in tumor-bearing mice. Activation of the p53 pathway-inhibited cell proliferation, as observed by detecting Ki67-positive cells. Additionally, the activity of apoptotic caspase family proteins was induced in the cleaved caspase-3-positive cells in vivo. Our results demonstrate that the MDM2 small-molecule inhibitor RG7388 is effective for NPC tumors, supporting further clinical investigation as a potential therapy for NPC.

Defective sarcomere assembly in smyd1a and smyd1b zebrafish mutants.

Two smyd1 paralogues, smyd1a and smyd1b, have been identified in zebrafish. Although Smyd1b function has been reported in fast muscle, its function in slow muscle and the function of Smyd1a, in general, are uncertain. In this study, we generated 2 smyd1a mutant alleles and analyzed the muscle defects in smyd1a and smyd1b single and double mutants in zebrafish. We demonstrated that knockout of smyd1a alone had no visible effect on muscle development and fish survival. This was in contrast to the smyd1b mutant, which exhibited skeletal and cardiac muscle defects, leading to early embryonic lethality. The smyd1a and smyd1b double mutants, however, showed a stronger muscle defect compared with smyd1a or smyd1b mutation alone, namely, the complete disruption of sarcomere organization in slow and fast muscles. Immunostaining revealed that smyd1a; smyd1b double mutations had no effect on myosin gene expression but resulted in a dramatic reduction of myosin protein levels in muscle cells of zebrafish embryos. This was accompanied by the up-regulation of hsp40 and hsp90-α1 gene expression. Together, our studies indicate that both Smyd1a and Smyd1b partake in slow and fast muscle development although Smyd1b plays a dominant role compared with Smyd1a.-Cai, M., Han, L., Liu, L., He, F., Chu, W., Zhang, J., Tian, Z., Du, S. Defective sarcomere assembly in smyd1a and smyd1b zebrafish mutants.

Gene Expression Program Underlying Tail Resorption During Thyroid Hormone-Dependent Metamorphosis of the Ornamented Pygmy Frog Microhyla fissipes.

Thyroid hormone (T3) is essential for vertebrate development, especially during the so-called postembryonic development, a period around birth in mammals when plasma T3 level peaks and many organs mature into their adult form. Compared to embryogenesis, postembryonic development is poorly studied in mammals largely because of the difficulty to manipulate the uterus-enclosed embryos and neonates. Amphibian metamorphosis is independent of maternal influence and can be easily manipulated for molecular and genetic studies, making it a valuable model to study postembryonic development in vertebrates. Studies on amphibian metamorphosis have been largely focused on the two highly related species Xenopus laevis and Xenopus tropicalis. However, adult X. laevis and X. tropicalis animals remain aquatic. This makes important to study metamorphosis in a species in which postmetamorphic frogs live on land. In this regard, the anuran Microhyla fissipes represents an alternative model for developmental and genetic studies. Here we have made use of the advances in sequencing technologies to investigate the gene expression profiles underlying the tail resorption program during metamorphosis in M. fissipes. We first used single molecule real-time sequencing to obtain 67, 939 expressed transcripts in M. fissipes. We next identified 4,555 differentially expressed transcripts during tail resorption by using Illumina sequencing on RNA samples from tails at different metamorphic stages. Bioinformatics analyses revealed that 11 up-regulated KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways and 88 Gene Ontology (GO) terms as well as 21 down-regulated KEGG pathways and 499 GO terms were associated with tail resorption. Our findings suggest that tail resorption in M. fissipes and X. laevis shares many programs. Future investigations on function and regulation of these genes and pathways should help to reveal the mechanisms governing amphibian tail resorption and adaptive evolution from aquatic to terrestrial life. Furthermore, analysis of the M. fissipes model, especially, on the changes in other organs associated with the transition from aquatic to terrestrial living, should help to reveal important mechanistic insights governing mammalian postembryonic developments.

Identification and differential regulation of microRNAs during thyroid hormone-dependent metamorphosis in Microhyla fissipes.

BACKGROUND: Anuran metamorphosis, which is obligatorily initiated and sustained by thyroid hormone (TH), is a dramatic example of extensive morphological, biochemical and cellular changes occurring during post-embryonic development. Thus, it provides an ideal model to understand the actions of the hormone and molecular mechanisms underlying these developmental and apoptotic processes. In addition to transcriptional factors, microRNAs (miRNAs) play key roles in diverse biological processes via post-transcriptional repression of mRNAs. However, the possible role of miRNAs in anuran metamorphosis is not well understood. Screening and identification of TH-responding miRNAs are required to reveal the integrated regulatory mechanisms of TH during metamorphosis. Given the specific role of TRs during M. fissipes metamorphosis and the characteristics of M. fissipes as an ideal model, Illumina sequencing technology was employed to get a full scope of miRNA in M. fissipes metamorphosis treated by T3. RESULTS: Morphological and histological analysis revealed that 24 h T3 treatment M. fissipes tadpoles resembled that at the climax of natural metamorphosis. Thus, small RNA libraries were constructed from control and 24 h T3 treatment groups. A total of 164 conserved miRNAs and 36 predicted novel miRNAs were characterized. Furthermore, 5' first and ninth nucleotides of miRNAs were significantly enriched in U in our study. In all, 21 miRNAs were differentially expressed between the T3 and control groups (p < 0.01). A total of 10,206 unigenes were identified as target genes of these differentially expressed miRNAs. KEGG pathway analysis indicated that the most overrepresented miRNA target genes were enriched in the "PI3k-Akt signaling pathway". In addition, a network associated with the TH signaling pathway provides an opportunity to further understand the complex biological processes that occur in metamorphosis. CONCLUSIONS: We identified a large number of miRNAs during M. fissipes metamorphosis, and 21 of them were differentially expressed in the two groups that represented two different metamorphic stages. These miRNAs may play important roles during metamorphosis. The study gives us clues for further studies of the mechanisms of anuran metamorphosis and provides a model to study the mechanism of TH-affected biological processes in humans.

Transcriptomics reveals the molecular processes of light-induced rapid darkening of the non-obligate cave dweller Oreolalax rhodostigmatus (Megophryidae, Anura) and their genetic basis of pigmentation strategy.

BACKGROUND: Vertebrates use different pigmentation strategies to adapt to various environments. A large amount of research has been done on disclosing the mechanisms of pigmentation strategies in vertebrates either under light, or, living in constant darkness. However, less attention has been paid to non-obligate, darkness dwellers. Red-spotted toothed toads Oreolalax rhodostigmatus (Megophryidae; Anura) from the karst mountainous region of southwestern China are non-obligate cave dwellers. Most tadpoles of the species possess transparent skin as they inhabit the dark karst caves. But remarkably, the transparent tadpoles can darken just within 15 h once exposed to light. Obviously, it is very significant to reveal molecular mechanisms of the unexpected rapid-darkening phenomenon. RESULTS: We compared the transcriptomes of O. rhodostigmatus tadpoles with different durations of light exposure to investigate the cellular processes and potential regulation signals for their light-induced rapid darkening. Genes involved in melanogenesis (i.e. TYR, TYRP1 and DCT) and melanocyte proliferation, as well as their transcriptional factor (MITF), showed light-induced transcription, suggesting a dominating role of morphological color change (MCC) in this process. Transcription of genes related to growth factor, MAPK and PI3K-Akt pathways increased with time of light exposure, suggesting that light could induce significant growth signal, which might facilitate the rapid skin darkening. Most importantly, an in-frame deletion of four residues was identified in O. rhodostigmatus melanocortin-1 receptor (MC1R), a critical receptor in MCC. This deletion results in a more negatively charged ligand pocket with three stereo-tandem aspartate residues. Such structural changes likely decrease the constitutive activity of MC1R, but increase its ligands-dependent activity, thus coordinating pigment regression and rapid melanogenesis in the dark and light, respectively. CONCLUSION: Our study suggested that rapid MCC was responsible for the light-induced rapid darkening of O. rhodostigmatus tadpoles. Genetic mutations of MC1R in them could explain how these non-obligate cave dwellers coordinate pigment regression and robust melanogenesis in darkness and light, respectively. To our knowledge, this is the first study that reports the association between pigmentation phenotype adaptation and MC1R mutations in amphibians and/or in non-obligate cave dwellers.

The complete mitochondrial genome of the Kaloula verrucosa (Anura: Microhylidae) and phylogenetic analyses.

The first complete mitochondrial genome of a Kaloula verrucosa frog was characterized in this work. The mitogenome was 17,061 base pairs (bp) in length, containing 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes, and a control region (D-loop). The overall base composition was 29.65% A, 30.77% T, 25.41% C and 14.17% G. Besides, the gene arrangement was identical to that observed in vertebrates. Five of 13 PCGs (COII, ATP6, COIII, ND3 and ND4) were ended with incomplete stop codon T. Except for ND6 gene encoded on L-strand, all other PCGs were encoded on H-strand. The non-coding region was 1665 bp in size, which was heavily biased to A + T (65.77%). Additionally, we found mitogenome size of all sequenced Kaloula species were bigger than that of Microhyla species, which was ascribe to the difference of D-loop size. Phylogenetic analysis showed that K. verrucosa was the sister species of Kaloula regifera. This work will provide basic molecular data for further molecular evolution and phylogenetic research of K. verrucosa and other microhylids.

De novo transcriptome assembly for the lung of the ornamented pygmy frog (Microhyla fissipes).

Microhyla fissipes, belonging to Neobatrachia, is a new model organism to study developmental biology, adaptive mechanisms from aquatic to terrestrial life, environmental toxicology, and human disease. M. fissipes use of lungs soon after hatching makes it extremely valuable for the study of lung function and development mechanisms. However, our knowledge of genes and pathways associated with lung development in M. fissipes is very limited. In this study, we conducted de novo transcriptome assembly for the lung of M. fissipes using the Illumina HiSeq4000 platform. We obtained approximately 9.0 GB clean data from the lung of the stage 28 tadpole with lung inflation (NCBI accession numbers: SRP107055). De novo transcriptome assembly identified 209,358 transcripts and 93,813 unigenes. In addition, BLASTX against NR, NT, KO, SwissProt, PFAM, GO and KOG databases were used to annotate all the 93,813 unigenes. This study provides the transcriptome and functional annotation of genes in M. fissipes lung development, which will be useful for comparative transcriptome analyses and promote research into mechanism of lung development in anuran.

Research proceedings on amphibian model organisms.

Model organisms have long been important in biology and medicine due to their specific characteristics. Amphibians, especially Xenopus, play key roles in answering fundamental questions on developmental biology, regeneration, genetics, and toxicology due to their large and abundant eggs, as well as their versatile embryos, which can be readily manipulated and developed in vivo. Furthermore, amphibians have also proven to be of considerable benefit in human disease research due to their conserved cellular developmental and genomic organization. This review gives a brief introduction on the progress and limitations of these animal models in biology and human disease research, and discusses the potential and challenge of Microhyla fissipes as a new model organism.

Transcriptome profiles of metamorphosis in the ornamented pygmy frog Microhyla fissipes clarify the functions of thyroid hormone receptors in metamorphosis.

Anuran metamorphosis is an excellent system in which to study postembryonic development. Studies on Xenopus (Mesobatrachia) show that thyroid hormone receptors (TRs) regulate metamorphosis in a ligand-dependent manner by coordinating the action of hundreds of genes. However, whether this mechanism is conserved among amphibians is still unknown. To understand the molecular mechanism of this universal phenomenon, we report the transcriptional profiles of the three key developmental stages in Microhyla fissipes (Neobatrachia): premetamorphosis (PM), metamorphic climax (MC) and completion of metamorphosis (CM). In total, 2,293 differentially expressed genes were identified from comparisons of transcriptomes, and these genes showed stage-specific expression patterns. Unexpectedly, we found that TRα was highly expressed in Xenopus laevis and Bufo gargarizans at premetamorphosis but showed low expression in M. fissipes. In contrast, TRß was highly expressed during metamorphosis in M. fissipes and X. laevis. This result may imply that TRß is essential for initiating metamorphosis, at least in M. fissipes. Thus, our work not only identifies genes that are likely to be involved in Neobatrachia metamorphosis but also clarifies the roles of unliganded TRα in regulating tadpole growth and timing of metamorphosis, which may be conserved in anurans, and the role of liganded TRß in launching metamorphosis.

The complete mitochondrial genome of Microhyla butleri (Amphibia, Anura, Microhylidae).

The complete mitochondrial genome was determined from a Microhyla butleri, Microhylidae, Microhyla, which was collected from Shenzhen, China. The mitogenome was 16 714 bp in length, containing 13 protein-coding genes, two rRNA genes, 22 tRNA genes and a control region (D-loop). The base composition was 28.7% A, 29.5% T, 27.2% C and 14.6% G. The gene order and contents were identical to most amphibian mitogenome. Except ND1 gene beginning with GTG and COI gene beginning with ATA, all other protein-coding genes began with ATG as start codon. Six protein-coding genes (ND1, COII, ATP6, COIII, ND3 and ND4) ended with incomplete stop codon T. The 22 tRNA genes with the size ranging from 65 bp to 74 bp were interspersed along the whole genome. The D-loop region containing tandem repetition was 1334 bp in length and heavily biased to A + T nucleotides.